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ABSTRACT 

A theory for the heating of coronal magnetic flux ropes is developed. The dissipated magnetic energy has 
two distinct contributions: (1) energy injected into the corona as a result of granule-scale, random footpoint 
motions, and (2) energy from the large-scale, nonpotential magnetic field of the flux rope. The second type of 
dissipation can be described in term of hyperdiffusion, a type of magnetic diffusion in which the helicity of the 
mean magnetic field is conserved. The associated heating rate depends on the gradient of the torsion parameter 
of the mean magnetic field. A simple model of an active region containing a coronal flux rope is constructed. 
We find that the temperature and density on the axis of the flux rope are lower than in the local surroundings, 
consistent with observations of coronal cavities. The model requires that the magnetic field in the flux rope is 
stochastic in nature, with a perpendicular length scale of the magnetic fluctuations of order 1000 km. 
Subject headings: MHD — magnetic fields — turbulence — Sun: corona — Sun: magnetic fields 
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1. INTRODUCTION 

In this paper we consider the heating of the solar corona 
in regions where the large-scale magnetic field deviates sig- 
nificantly from a (current-free) potential field. Evidence for 
such nonpotential structures comes from a variety of sources. 
Measurements of photospheric vector fields in act ive regions 
often show deviations from the potential field (e.g ., lGarv et alj 
fl987t lPevtsovetal][T99^. fl997t iLeka et alll996h . X-ray ob- 
servations show sigmoids, which are bright S-shaped or in- 
verse S^shjvpe^lstructOT active regions 
(e.g.. lActon et all 1 19921: iRust & Kumarlll996l ICanfield et alj 
Il999t ISterling"e t al. 2000; Gibson j |20025~These sigmoids 
contain coronal loops that are highly sheared with respect to 
the photospheric polarity inversion li ne (PIL), and are are of- 
ten associated with Ha filaments (e. g., Manoharan et a ll 1996) : 
lPevtsovl200"llGibson & Fanl2006a ). Schmi eder et all i ll 9961) 
proposed a differential magnetic shear model in which the 
degree of s hear decreases with h eight in the observed ac- 
tive region. Schrij ver et alJ ((2005) studied the nonpotential- 
ity of active regions using extreme-ultraviolet (EUV) images 
from the Transition Region and Coronal Explorer (TRACE). 
They found that significant nonpotentiality occurs when new 
magnetic flux has recently emerged into the corona, or when 
rapidly evolvin g, opposite polarity flux concentrations are in 
close contact. Mand rini et alj ([2000) tested various models of 
coronal heating using potential, linear force-free, and magne- 
tostatic models. Of course, nonpotential magnetic fields also 
play an i mportant role in sola r flares and coronal mass ejec- 
tions (see lFan & Gibsor]|2007l and references therein). 

The nonpotential structures in the corona have been mod- 
eled as flux ropes in which the field lines twist about an ax 
ial field line (e. g., Rust & Kumarlll994t lLow & H undhause 
19951: lAulanier & Demoulinlll998l: iTitov & Demoulinlll99 



Amari et al. 2003; Magara & Longcope 2001; Kliem et al 



Kusano 2005; Aulan ier et all 120051; iGibson et al 



The presence of a coronal flux rope implies the 
existence of electric currents that pass through the photo- 
sphere and into the corona. The corona of an active re- 
gion is a magnetically dominated plasma. Therefore, ex- 
cept during a flare or eruption, the corona is nearly force 
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free, VxBk aB, where B(r) is the magnetic field and a(r) 
is the so-called torsion parameter. A magnetic field with 
this property is called a nonlinear force-free field (NLFFF). 
Many authors have developed NLFFF models of active re- 
gions, either by extrapolating photospheric vector fields into 
the corona (e . g.. IWheatland et alJl2000t iBlevbel et al l [2002 : 
Regnier et all 120021: IWiegelmannf 120041; IValori et all 120051 
Wheat land 2006; Regnier & Priest 2007), or using a flux rop e 
insertion method (van Ballegooiien 20041 iBobra et alj 2008). 

Active regions exhibit transient brightenings at EUV and 
soft X-ray wavelengths on time scales of a few minutes 
(e.g.. IShimizu et all [l 992h iShimizul 1 1 995b iPorter et al.lll995b 
Berghman s"& Clettel Il999t iBerghmans et all 120011) . These 
single- or multi-loop brightenings are thought to be due to 
the release of magnetic energy by small reconnection events 
in the complex coronal field of an active region. This may be 
a two-step process: magnetic energy may first be converted 
into kinetic energy of reconnection outflows, waves or en- 
ergetic particle s, and then be c onverted into thermal energy 
of the plasma (Longcope 2004). It seems likely that X-ray 
sigmoids are also heated by such reconnection events, with 
many events occuring simultaneously at different sites within 
the sigmoid. The events may occur so frequently that the 
coronal plasma does not have time to cool significantly be- 
tween heating events, keeping the plasma at temperatures of 
3-5 M K (e.g.. lKlimchuk & Porter|[l99l ICargill & Klimchukl 
120041) . This would explain why sigmoids generally do not 
show up in EUV spectral lines formed at 1-1.5 MK. 

What is the origin of the energy released in small reconnec- 
tion events? One possibility is that the energy is injected into 
the corona by small-scale random footpoint motions, which 
lead to twist ing and braiding of coronal field lines (e.g- lParkerl 
1 19721 1 1988b . However, in the highly sheared magnetic field 
of a sigmoid there is another possibility, namely, that the en- 
ergy for coronal heating is derived from the energy of the flux 
rope itself. In this case the energy is already stored in the 
corona and does not need to be injected at the photospheric 
footpoints. Generally, a large amount of free energy is avail- 
able in the coronal flux rope (~ 10 32 erg), sufficient to explain 
both the heating of the coronal plasma and the occurrence of 
flares and coronal mass ejections. 

The second scenario leads us to propose a new theory for 
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FIG. 1 . — Stochastic magnetic fields in a coronal flux rope. 



the heating of coronal flux ropes. The basic ideas are illus- 
trated in Figure Q] The lower part of the figure shows a flux 
rope in the low corona. Inside the flux rope is a magnetic 
field that runs more or less parallel to the flux rope axis. The 
flux rope is held down by an overlying coronal arcade (not 
shown), which has a magnetic field that runs perpendicular to 
the flux rope axis. Hence, going from the inside to the out- 
side of the flux rope the direction of the mean magnetic field 
rotates through an angle of about 90 degrees. The upper part 
of the figure shows the fine-scale structure of the magnetic 
field at the boundary between the flux rope and its surround- 
ings. The magnetic field has both mean and fluctuating com- 
ponents, as indicated by the wiggles in the field lines. The 
total magnetic field B(r) is as sumed to be stochastic in na- 
ture (lLazarian & Vishniaclll999l) . The magnetic fluctuations 
are assumed to be produced by some unspecified MHD pro- 
cess that ultimately derives its e nergy from the free magnetic 
energy of the coronal flux rope. Strauss! (119881) proposed that 
turbulent fluctuations could be produced by tearing modes in 
coronal current sheets, but the present paper does not depend 
on that particular mechanism. We only assume that reconnec- 
tion occurs intermittently at many sites in the stochastic field. 
The reconnection events produce outflows and energetic par- 
ticles that ultimately result in heating of the coronal plasma. 

The above process is an example of a "dynamo" pro- 
cess because it involves turbulent diffusion of the mean 
magnetic field of the coronal flux rope. Dynamo ac- 
tivity, defined as the amplification and/or maintenance of 
large-scale magnetic fields, occurs in many astrophysical 
and laboratory plasmas, and magnetic he licity evolution 
plays an important role in such activity (e .g., Blackm an & Jil 
2006; Subrama rrian & Brandenburg! ]2006l) . iBoozerl (119861) 
and iBhattachariee & Hameiril d!986l) considered turbulence 
in low-/? plasmas, and proposed that the evolution of the 



mean magnetic field can be described in term of hyperdif- 
fusion, a type of magnetic diffusion in which the ma gnetic 
helicity of the mean field is conserved. IStraussI d!988l) pro- 
posed that tearing modes can produce turbulence in three- 
dimensional sheared magnetic fields, and showed that such 
turbulence causes hyperdiffusion, which may explain both 
fast recon nection in solar flares and th e heating of the so- 
lar corona. Bhattach ariee & Yuanl (Tl995) considered the con- 
straints on the dynamo mechanism from helicity conserva- 
tion and the back-reaction of the magnetic field on plasma 
dynamics. They showed that the results of kinematic dy- 
namo theory are strongly modified by the production of hy- 
perdiffusion in the mean field equation. Hyperdiffusion has 
also been included in models for decaying active regions 
( Ivan Ballegooiien & Mackavll2007l) . 

In the present paper we consider the plasma heating asso- 
ciated with hyperdiffusion. At present, it is unclear to what 
extent stochastic effects dominate the solar magnetic field. 
Our purpose here is not to develop a complete theory of such 
fields, but rather to determine some of the observable con- 
sequences of their existence. Therefore, more detailed ques- 
tions, such as how the stochastic field is produced and the 
nature of the MHD instabilities, will remain unanswered. The 
main goal is to demonstrate that the existence of stochastic 
fields is consistent with existing observations of coronal heat- 
ing in active regions. 

2. THEORY OF HYPERDIFFUSION 

The theo ry of hyperdiffusion was f i rst pro posed by lBooze"rl 
(1986) and Bhattach ariee & Hameiril d 19861). and fur t her de - 
veloped by IStraussI d!98 8). Bha ttachariee & Yuan! d 19951) . 
Vishni ac & Choi (1200 ll) and others. In this section we sum- 
marize the results of these analyses as they apply to the solar 
corona, i.e., when the plasma is magnetically dominated. The 
total magnetic field B(r,f) is a sum of mean and fluctuating 
components, and is described in terms of the vector potential 
A(r,f), where Be VxA. The evolution equation for the 
vector potential is 

BA 

— =-cE-V$, (1) 
at 

where E(r,f) is the electric field, Q(r,t) is the scalar potential, 
and c is the speed of light. Inserting Ohm's law and Ampere's 
law, we obtain 

<9A 

— =vxB-nVxB-V$, (2) 

at 

where v(r,f) is the plasma velocity, and rj is the Ohmic diffu- 
sivity. We assume that ij <C L 2 /t, where L is the system size 
and r is the time scale for the evolution of the mean magnetic 
field. Then reconnection occurs only in highly localized re- 
gions (e.g., current sheets ) that make up only a sm all fraction 
of the volume. Following I Vishniac & Choi ([200 1 ), we use the 
Coulomb gauge, V • A = 0, so that the scalar potential is given 
by 

V 2 $ = V-(vxB-ijVxB). (3) 
The magnetic helicity equation follows from equation (|2]): 

4-(A-B)+V-W = -2?7B- VxB, (4) 

at 

where (A • B) is the helicity density, and 7i(r, t ) is the helicity 
flux: 

HeAx(vxB-t)VxB+V$) 

= Ax(-cE+V$). (5) 
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Bergerj d 19841) has shown that in the limit rj <C L 2 /t the right- 
hand side of equation can be neglected, so this equation 
describes the conservation of magnetic helicity. The relative 
magnetic helicity is defined by H« = /[A-B- A p - Bp] dV, 
where subscript p refers to the potential field, and the in- 
tegration includes th e coronal and subsurface regions (see 
lBerger&Fieldlll984l) . Therefore, in the limit 77 < L 2 /t, if 
no helicity is injected into the corona at the photosphere or 
ejected into the heliosphere, Hr is constant in time. 

Let Bo(r, t ) be the mean magnetic field smoothed over some 
intermediate spatial scale i, and let Bi(r,f) be the fluctuating 
field on scales less than £. Similar decompositions are made 
for other physical quantities. Applying the smoothing opera- 
tor to equations (Q~|i and (0, we obtain: 

<9A 

— - = -cE - V$ « y'o x B + £ - V* , (6) 
at 

where Ao, &o, Eo and v are the mean vector potential, 
scalar potential, electric field and plasma velocity, and £q=< 
Vi x Bj > is the "mean electromotive force" describing 
the effect of the flu ctuations on the mean magnetic field 
dSteenbeck et al.l [19661) . Here we neglected the Ohmic dif- 
fusion of the mean field. The helicity equation for the mean 
magnetic field can be derived from equation (0: 

^ (Ao • B )+ V • [Ao x (t x Bo + £ Q + V* )] = 2fb ■ Bo . (7) 

On the other hand, applying the smoothing operator directly 
to equation with 77 « 0, we obtain: 

-|«A-B» + V- [A x (v xB +So+V$o)+Wo] =0, 

(8) 

where Ho is the contribution of the fluctuations to the helicity 
flux: 

H = <Ai x(-cE 1 + V$ 1 )> . (9) 

We assume that the large-scale structures dominate the he- 
licity budget, so that the helicity of the total field can 
be approximate d by the he l icity of the mean field, < A • 
B A ■ Bp (Boozer| [T986HBhattachariee & Hameirill 19861: 
IVishniac & Choll2001l) . Subtracting equation (7} from equa- 
tion ([8]l yields the following relationship between £0 and Ho: 

2£o-B = -V-H . (10) 

The electromotive force can be written as a sum of compo- 
nents parallel and perpendicular to the mean magnetic field, 
£0 = £pj| +£o,_i_- The parallel component is given by equa- 
tion ( 1101 ). and the perpendicular component can be written in 
terms of a transport velocity uo, which is defined by uo = 
(Bo x £o j_)/Bq. Then the electromotive force is given by 

£o = u x B - \ J^° Bo. (11) 

We assume that the magnetic field evolves slowly in re- 
sponse to hyperdiffusion and changes in the boundary condi- 
tions. Furthermore, we assume that the magnetic pressure is 
much larger than the plasma pressure, so that the mean mag- 
netic field evolves through a series of force-free equilibrium 
states, (V x Bo) x Bo « 0. The force-free condition can also 
be written as 

VxB (1 = q„B , (12) 

where ao is the torsion parameter of the mean field. For a lin- 
ear force free field (ao = constant) the helicity flux Ho should 



vanish because there is no magnetic free energy avail able to 
drive t he m agnetic fluctuations. Ther e fore, f ollowing iBoozerl 
( 1986) and Bhattacharjee & Hameiri ( 1986) we assume that 
Ho is proportional to the gradient of ao: 

Ho=-2 m B 2 Va , (13) 

where 774 is the so-called hyperdiffusivity. Inserting expres- 
sion dT3l into equations ( fTTT i and (0 yields the following dif- 
fusion equation for the mean magnetic field: 

= v X B + ^ V ■ ( m BlVao)- V$ , (14) 

where Vo = V + Uo is the total transport velocity. The sec- 
ond term on the right-hand side of equation dT4b describes the 
effect of the magnetic fluctuations in diffusing the mean mag- 
netic field. The effect is such that the magnetic helicity of the 
system is conserved. Hyperdiffusion describes the t endency 
of the magnetic field to relax to a state of constant ao dTavlor I 
119741) . The heating of the coronal p lasma resulting from Tay - 
lor relaxation was first modeled bv lHevvaerts & Priest! d!984l) . 

The energy equation for the mean magnetic field can be 
derived from equations (O, dTTT l and dT2i >: 

(15) 

and inserting equation ( fT3T > we can write 

s(S) +v '(S^-*S a »S=-*S !V<w|2 

(16) 

where vo,j_ = v j_ +Uo is the total transport velocity. The quan- 
tity on the RHS of equation dTST ) is negative and can be inter- 
preted as t he rate at which mean magnetic energy is con verted 
into heat dBoozed l!986: Bhatt acharjee & Hameirilll986l) . As- 
suming there is no inflow of energy at the boundaries, the total 
magnetic energy decreases monotonically with time: 

3. HELICITY FLUX IN STOCHASTIC MAGNETIC FIELDS 

The above theory predicts a particular form of the mean- 
field diffusion equation, but does not explain why the helicity 
flux Ho is proportional to Vao, nor does it predict the value 
of the hyperdiffusivity 774(1", f) contained in equation dT3l l. To 
estimate this quantity, a detailed understanding of the pro- 
cesses that driv e the magnetic fluctuations would be required. 
IStraussi d!988l) assumed the fluctuations are driven by tear- 
ing mode instability, but in the present paper we make no as- 
sumption about the nature of the instability. We only assume 
that ( 1 ) the driving process results in stochastic magnetic field 
lines dLazarian & Vishniacll 1999b . as illustrated in Fig. [TJ; (2) 
reconnection occurs intermittently at many sites within this 
stochastic field; (3) the reconnection events result in Alfvenic 
perturbations that travel along the stochastic field lines until 
they are damped as a result of turbulent cascade. Therefore, 
the magnetic field is not completely force free. In a non-force- 
free field, the magnetic torsion parameter a can be defined as 
a = (B • V x B)/B 2 , and we assume that this quantity varies 
in a time-dependent way along the field lines. The topology 
of the field is continually changing as a result of the reconnec- 
tion processes, with new magnetic connections being formed 
and other connections being severed. 
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Consider a small volume near the edge of the coronal flux 
rope where the mean torsion ao(r) can be approximated as a 
linear function of position. Let z be the coordinate perpendic- 
ular to the surfaces of constant ao (we assume dao/dz > 0). 
We assume that at time t = magnetic reconnection occurs at 
some localized site, creating a new magnetic connection be- 
tween two neigboring surfaces z = constant. Let Pi and P 2 
be two points on the newly formed field line, such that P\ is 
located on z = Z\ and P 2 is located on z = Z2 > Z\- The dis- 
placement As between these two points as measured along 
the field line is assumed to be much larger than the distance 
Az between the two surfaces. The displacements are given by 



Az ~ Aj 



A*~ A 



I- 



(18) 



where An and X± are the parallel- and perpendicular length 
scales of the turbulent motions (relative to the mean mag- 
netic field). The torsion parameters a(P\) and a(P2) at the 
two points differ by an amount 



Aa = a(P 2 )-a(Pi) « aofo)-ao(zi) ~ Az 



dz 



(19) 



We assume that this difference in torsion results in the launch 
of Alfven waves that travel in both directions along the newly 
formed field line. In Appendix A we show that these waves 
carry a certain amount of magnetic helicity. The net helicity 
flux due to the waves can be estimated by inserting expression 
( IT9l > into equation (lA14b : 



H^-f \\v A Bl Az^ 
dz 



(20) 



where v A is the Alfven speed, and we assume that the radius R 
of the newly formed flux tube is determined by the outer scale 
of the turbulence, Rk X±. The quantity /q is a dimensionless 
constant, and we estimate it to be rather small, /o « 0.01 (see 
Appendix). The z-component of the helicity flux is 

Az 
As 



m^K-foxWv A Bl^ 



dz 



(21) 



where we take into account the angle of the field lines with 
respect to the z-direction, and we use equations < f20b and ( fT8l . 
This shows that the helicity flux Ho due to the fluctuations is 
proportional to the local gradient of ao, as assumed in equa- 
tion ( fT3] >. Furthermore, we obtain the following estimate for 
the hyperdiffusivity: 



r/4~ foX 4 ± X^v A . 



(22) 



Note that 774 depends strongly on the perpendicular length 
scale of the turbulent motions associated with the reconnec- 
tion in the stochastic magnetic field. 

4. CORONAL HEATING 

The coronal magnetic field is anchore d in the photo sphere 
where the field is highly fragmented dStenflol fl989). The 
photospheric field consists of vertically oriented, kilo-Gauss 
flux tubes separated by nearly field-free plasma. These flux 
elements are continually buffeted by convective flows on 
the scale of the solar granulation, which results in random 
twisting of the flux tubes and a rand om walk of the flux 
elements over the so lar surface (e.g., [Berger & Title! 1 19961; 
iNisenson et al.l p2003). Both types of motions cause distur- 
bances that propagate upward along the field lines in the f orm 
of Alfven waves (e.g. JCranmer & v an Ballegooii enl2005h . In 
coronal loops, such disturbances produce small-scale twisting 



and braiding of the coro nal field lines and associated fi eld- 
aligned electric currents ( Parker] [19721) . ISchriiverl (120071) ar- 
gues that the observed widths of coronal loops are determined 
by these footpoint motions. The random walk of the foot- 
points is predicted to drive a ca scade of magnetic energy from 
large to small spatial scale s (Ivan Ballegooiien 1985, 11986L 
119881: ICraig&Snevdl 120051) . The existence of such a cas- 
cade is confirmed by results from 3D numerical simulations of 
the Parker problem with incompressible boundary flows (e.g., 
i Mikic et al.lll989HHendrix & van Hpv"enlll996l:lHendrix et al.l 
1 19961: iGalsgaard &Nordlundl 119961: iRappazzo et alj 120071) . 
However, the twisting motions of the flux tubes may be more 
important than the lateral motions, as the twi st velocity is am- 
plified by the expansion of the flux tubes (Ivan Ballegooiienl 
1986; Man drini etaO 120001) . At the smallest scales in the 
corona, impulsive reconnection events ("nanoflares") play 
an important role in the release o f the magnetic energy 
dParkerl 119881; iLu & Hamilton! 1199 lb . The actual dissipa- 
tion may be a multi-step process invol ving reconnectio n 
outflows, waves, and energe tic particles (Long copel 120041) . 
iGudiksen & Nordlundl (120051) have modeled the coronal heat- 
ing that ultimately results from granule-scale footpoint mo- 
tions. 

In this paper we assume that the coronal heating rate e has 
contributions from two different sources: (1) disspation of 
energy that was directly injected into the corona by granule- 
scale footpoint motions, as described above; (2) dissipation of 
energy from a large-scale coronal flux rope via the process of 
hyperdiffusion. The former uses energy propagated upward 
along the field lines, while the latter draws its energy from 
large-scale magnetic shear already present in the corona. The 
direct contribution to the heating is approximately given by 



. Bf^u^s)^ 



(23) 



where s is the position along the coronal loop, edi rec t(s) is 
the heating rate per unit volume, Bq(s) is the magnetic field 
strength, u(s) is the velocity of the perpendicular motions, r p h 
is the correlation time of the photospheric motions, L is the 
loop length, and < ? » 3. This expression is a generalization 
of equation (61) in lvan Ba llegooiien ( 1986), using a magnetic 
Reynolds number R m ~ 10 ll) . Equation ( 1231 assumes that the 
Alfven-travel time along the coronal loop is short compared to 
the correlation time, L/v A < r p h. The perpendicular velocity 
is given by 



U(s) f» H p h 



B 



P h 



B (s)' 



(24) 



where M p h is the root-mean-square velocity in the photosphere, 
Bph is the photospheric field strength, and we assume that the 
magnetic flux associated with the perturbations is constant 
along the field lines from the base of the photosphere into the 
corona. Combining the above equations, the direct heating 
rate is 

. B ph«ph r ph 



£direct(i) « q 



4irL 2 



-Bo(s), 



(25) 



similar to equation (67) of Ivan B allegooiien (1986). In 
the present paper we use B p h = 1 kG, u p h = 0.35 km s -1 , 
and r p h = 600 s, which yields a photospheric diffusion con- 
stant D = 92 km 2 s" 1 [see equation (62) of Ivan B allegooiien 
( 1986)], consistent with the observed dispersal coefficient that 
characterizes the gran ular random walk up to several hours 
dHagenaaretal.ll 1997b . 
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The heating rate due to hyperdiffusion follows from equa- 
tions dT6j and d22l >: 

, B 2 (s)X 4 ± (s) 2 
'fo-^—T — v A (s)\Va \ , (26) 



EhyperCs) = 774^ | Va | 2 



4ttL 



where we assume that the parallel length scale of the fluctua- 
tions is determined by the loop length, Ay ~ L. The perpen- 
dicular length scale A j_Cs) of the turbulence is assumed to vary 
along the field lines in a manner consistent with magnetic flux 
conservation: 



B 



Ai(i)«W-f-, (27) 

where A tm b is the corresponding length scale in the photo- 
sphere. However, the turbulence is not driven by the foot- 
point motions, therefore A tur b does not necessarily reflect any 
property of the photospheric velocity field. We cannot pro- 
vide a more detailed expression for A tul -b because the nature 
of the physical processes that produce the stochastic field are 
still unclear; instead, we treat A tur b as a free parameter of the 
model. Combining the above expressions, we find 



Chyperfa) ~ /o 



R2 \4 
B ph A turb 

4ttL 



va(j)|Voo| : 



(28) 



The total heating rate is assumed to be the sum of expressions 
(ESl and (|28): 

e(s) = £diiect(«) + EhyperW- (29) 

In the next section we use these expressions to compute the 
coronal temperature and density, which provides a link to so- 
lar observations. 

5. ACTIVE REGION MODEL 

iBobra et ail (120081) developed empirical models for 
NLFFFs in two active regions observed with TRACE. Our ini- 
tial goal was to compute temperatures and densities for these 
models. However, we found that small deviations from the 
force-free condition can cause significant errors in the dis- 
tribution of Veto along the field lines, so that the heating 
rate ehy P er(s) cannot yet be accurately computed from such 
empirical models. Therefore, in the present paper we apply 
the above expressions for the heating rate to a semi-analytic 
model of an active region. The model describes the mean 
magnetic field Bo(r) in the corona; fluctuating fields are not 
modeled, nor do we consider the evolution of the mean field. 

The NLFFF model is constructed as follows. We use a 
Cartesian coordinate system (x,y, z) with z being the heigh t 
above the coronal base. Following lTitov & Demoulinl(ll999t) . 
we assume that the magnetic field has cylindrical symmetry 
with respect to an axis that lies horizontally below the coro- 
nal base at y = 0, z = —R. Let (r, <j>,x) be a cylindrical coor- 
dinate system with r the distance from the cylinder axis and 
x the distance along the axis. We consider a finite domain 
r = [R,R m ] and x = [— x m ,x m ], where the inner radius of the 
domain is given by R, and x m 3> R- We initially use dimen- 
sionless units such that R = 1 . The force free field is written 
as: 

B(A) 

(30) 



1 OA I OA 

Bo.r = 5-, #o.a = - -5- , 

r ox r or 



Bo,<p = 



where A(r,x) is the flux function, and B(A) describes the az- 
imuthal component of the field. The flux function satisfies the 
following non-linear equation: 



d ( 1 dA \ d 2 A 
i" tt ( ) + "o~T + a o(A)B(A) = 0, 
Or \ r or I ox z 



(31) 




FIG. 2. — Non-linear force-free field model of an active region contain- 
ing a magnetic flux rope. The contours show the flux distribution on the 
photosphere (red: positive flux; green: negative flux). The black lines show 
selected field lines in the flux rope. 

where ao(A) = cLB/dA. In this paper we use boundary con- 
ditions such that A(l,x) = 0.8exp(-x 2 ) at the inner boundary, 
and A = at the outer boundaries. The azimuthal field is given 
by 

B(A) = (3[l-exp(-A 5 )]. (32) 

where (3 is a parameter describing the magnitude of the az- 
imuthal field. Equation PIT ) is solved using a relaxation tech- 
nique. In essense, we solve the following time-dependent 
equation: 



dA 
~dt 



■ D 



d_ 
dr 



1 dA 
r dr 



d 2 A 

a? 



+ a (A)B(A) 



(33) 



until a stationary solution A(r,x) is obtained (here D is a di- 
mensionless diffusion constant). In each time step, (3 is ad- 
justed such that the total azimuthal current has a fixed value, 
/ / ccqB dxdz = 2.5. Once a stationary solution is obtained, 
the spatial scale is changed to R = 30 Mm, and the magnetic 
field is rescaled by a factor 100, producing a maximum value 
of \B()^\ of 55 G. The quantities Bo.r, #0,* and Bq^ are com- 
puted on a grid of 200 x 200 points in the {r,x) plane. The 
magnetic vector Bo(x,y,z) in the cartesian frame can then be 
found by evaluating these quantities at r = \fy 2 + {R+z) 2 . Fig- 
ure [2] shows the magnetic flux distribution at the coronal base 
(z = 0) and selected field l ines in the flux rope. The stability 
of the iTitov & Demoulinl i 19991) model has been considered 
bv lTorok et al.l (120041) . They find that the m = 1 kink mode is 
stable provided the average twist in the flux rope is less than 
3.57T. This condition is satisfied in the present case, so we 
assume that the field shown in Fig.[2]is stable. 

The NLFFF model provides the spatial distribution of 
cto(r), which is needed to compute the heating rate e. We 
trace a large number of field lines that cross the vertical plane 
y = 0, and solve the energy balance equation for each field line. 
W e use a modified version of the iterat ive method described 
by ISchrijver & vanB allegooijen (2005). The code assumes 
steady-state heating, and includes thermal conduction and op- 
tically thin radiative losses. The contribution to the heating 
from hyperdiffusion is determined from equation (f28b with 
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FIG. 3. — Results from modeling the energy balance of the plasma along 
a particular magnetic field line located at the edge of the coronal flux rope 
where | V«o I is largest, (a) Temperature T(s) (full curve) and electron density 
n e (s) (dashed curve) as functions of position s along the field line, (b) Total 
heating rate t(s) (full curve), which is a sum of hyperdiffusive heating and 
a "direct" contribution from photospheric footpoint motions. The latter is 
shown by the dashed curve. 



the Alfven speed Va{s) taken from the previous iteration. We 
assumed that the perpendicular length scale of the turbulence 
is characterized by A tur b = 10 3 km. The chromosphere-corona 
transition region (TR) is assumed to be located at z = 0, where 
the temperature T = 2 x 10 4 K. The solution of the energy bal- 
ance equation yields the temperature T{s) and electron density 
n e (s) as functions of position s along a field line. 

Figure [3] shows modeling results for one field line at the 
edge of the coronal flux rope. Starting at x = y = and z = 43 
Mm, the field line is traced forward and backward through the 




FIG. 4. — Coronal models for an active region containing a magnetic flux 
rope. The upper panels show results for a model with both direct and hyper- 
diffusive heating: (a) temperature T(x,z) in a vertical plane y = intersecting 
the flux rope; (b) electron density n e (x,z) in the same vertical plane. The 
height and width of these grey scale images correspond to 100 Mm on the 
Sun. The maximum temperature in the model is 2. 1 1 MK, and the maximum 
density is 2.7 X 10 s cm" 3 . The bottom panels show results for a model with- 
out hyperdiffusion: (c) temperature T(x,z), and (d) electron density n e (x,z)- 
This model has lower peak temperature and density (maximum values are 
1.42 MK and 1.0 X 10 9 cm" 3 , respectively). In both models the minimum 
temperature is 0.37 MK, and the minimum density is 1.8 X 10 7 cm" 3 . Maxi- 
mum values are shown as white; minimum values are black. 



3D magnetic model until the level z = is reached. The tem- 
perature T(s) and density n e (s) along this field line are shown 
in Fig. [3^, and the heating rate e(s) is shown in Fig. (3j> (full 
curve). The dashed curve in Fig. shows the "direct" con- 
tribution to the heating from photospheric footpoint motions. 
Note that this contribution is relatively small, so for this par- 
ticular field line the heating is dominated by hyperdiffusion. 
However, for other field lines in the NLFFF model (e.g., along 
the axis of the flux rope) the hyperdiffusive heating rate is 
smaller than the direct heating rate. 

We computed temperature and density profiles for a large 
number of field lines that cross the vertical plane y = 0. Fig- 
ure|4]shows the temperature T(x,z) and density n e (x,z) in the 
plane y = for two different heating models. The top panels 
show the results when both direct and hyperdiffusive heating 
are applied together (A tur b = 10 3 km). Note that the tempera- 
ture and density are enhanced in a thick shell where the gradi- 
ent of a.Q is largest. The density distribution shows a "cavity" 
that is round at the top and V-shaped at the bottom. The tem- 
perature is also enhanced at the bottom of the V-shaped cavity 
because the field lines in this region are relatively long, so the 
conductive cooling is reduced. 

The bottom panels of Figure [4] show the temperature and 
density distributions in the y = plane when only the direct 
heating is applied (i.e., no hyperdiffusion). The highest tem- 
peratures now occur within the flux rope where the field lines 
are longest and conductive losses are weakest. The lowest 
densities occur in the V-shaped lower part of the flux rope, 
but the density in the upper part is nearly the same as in the 
surroundings, and there is no dome-shaped cavity. 
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On the Sun, cavities are frequently observed in the K- 
corona emission, and these cavitie s have been linked to the 
presence of coronal flux ropes (see iGibson et al.ll2006bl and 
references therein). We suggest that the observed cavities are 
evidence for a spatial distribution of coronal heating that is 
consistent with hyperdiffusion. 

6. DISCUSSION AND CONCLUSIONS 

We propose a new mechanism for coronal heating relevant 
for magnetic structures in which the large-scale magnetic field 
is highly sheared and/or twisted. We assume that the mag- 
netic field in such coronal flux ropes is stochastic in nature 
dLazarian & Vis hniac 1999), and contains small-scale pertur- 
bations superposed on a mean magnetic field. The magnetic 
free energy associated with the flux rope can be converted into 
heat via small-scale reconnection processes. We describe this 
process in terms of hyperdiffusion, a type of magnetic diffu- 
sion in which the magnetic helicity is conserved. The effect 
of hyperdiffusion on the mean field depends on the gradient 
of the mean torsion parameter ao(r,f). Therefore, hyperdiffu- 
sion describ es the tendenc y of a magnetic system to relax to a 
Taylor state dTavlor II 1974b . 

In section 3 we consider helicity transport in a stochastic 
magnetic field, and we show that small-scale reconnection 
events produce a net helicity flux Ho consistent with hyper- 
diffusion. We derive an expression for the hyperdiffusivity 7/4 
in terms of the parallel and perpendicular length scales of the 
magnetic perturbations. The calculation is based on a sim- 
ple model for the Alfven waves produced by reconnection 
events (Appendix A). We find that the hyperdiffusivity de- 
pends strongly on the perpendicular length scale of the turbu- 
lence, j]4 oc AjL b . The spatially averaged heating rate eh ype r due 
to hyperdiffusion is proportional to r^lVaol 2 . This heating is 
distinct from the "direct" heating due to small-scale random 
footpoint motions. In section 4 we derive an expression for 
the total heating rate e, which is a sum of direct and hyperdif- 
fusive contributions. 

Our analysis of the dissipation process is very approximate 



and makes many simplifying assumptions. For example, we 
cannot identify a specific mechanism by which the stochas- 
tic magnetic field is produced and maintained. We simply 
treat the perpendicular length scale A tlu b of the fluctuating 
field as a free parameter to be determined from observations 
or more complete theories of turbulent reconnection. Also, 
it is unclear exactly what drives the small-scale reconnection 
events necessary for hyperdiffusion. The photospheric mag- 
netic field is highly fragmented, and there are separatrix sur- 
faces in the corona corresponding to the interfaces between 
photospheric flux tubes. We speculate that the reconnection 
events may occur at these separatrix surfaces. 

In section 5 we develop a 3D model of an active region 
containing a magnetic flux rope, and we compute the temper- 
ature and density distributions with and without hyperdiffu- 
sion. The results show that hyperdiffusion can have a signifi- 
cant effect on the temperature structure, provided the perpen- 
dicular length scale of the magnetic fluctuations A tllr b ~ 10 3 
km at the coronal base. If A tur b were of order 100 km or less, 
the hyperdiffusive heating would not have a significant impact 
on the total heating rate for any field line. For the NLFFF 
model considered here, the heating is concentrated in a thick 
shell surrounding the flux rope. The heating at the axis of the 
flux rope is weaker, producing a "cavity" similar to the coro- 
nal cavities observed on the Sun. 

At present there is no direct observational evidence to con- 
firm or refute the existence of stochastic magnetic fields with 
perpendicular scales ~ 10 3 km. Therefore, it is unclear ex- 
actly how coronal flux ropes are heated. Coronal images ob- 
tained with the Solar Dynamics Observatory (SDO) may pro- 
vide new observational constraints on stochastic fields and the 
physical processes by which coronal flux ropes are heated. 
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APPENDIX 

TRANSPORT OF HELICITY BY ALFVEN WAVES IN A STOCHASTIC FIELD 

The transport of magnetic helicity plays an important role in the theory of hyperdiffusion. Here we consider a simple model 
for helicity transport due to reconnection events within the stochastic field. The average time T between reconnection events on a 
given flux tube is assumed to be a few times the Alfven travel time along that tube. Let us assume that at time t = a new section 
of magnetic field has just been created by reconnection. The new section connects two points P\ and Po that may be located at 
opposite ends of a coronal loop. We neglect the loop curvature and approximate the new section as a cylindrical flux tube with 
radius R. We assume that in the initial state the plasma in the tube is stationary, but the torsion parameter a varies along the tube, 
so the field is not force free. The variable torsion causes the launch of counter-propagating Alfven waves that travel in opposite 
directions along the tube. We assume that the wave length Ay of the waves is large compared to the tube radius. The medium 
surrounding the flux tube is assumed to be unaffected by the waves. 

We use a cylindrical coordinate system (r, (fi,s) where s is the distance along the tube (0 < s < As). The plasma is assumed to 
be fixed at the end points Pi (s = 0) and P2 (s = As), so that the waves reflect at these points. The magnetic field B inside the tube 
is written as 

E(r,s,t)=B s + ^ [BfrQ + SB^t)] <p, (Al) 

where Bq is the axial field (constant in space and time), B^jo is the unperturbed azimuthal field at the cylinder wall (r = R), and 
SB^s^t) is the perturbation of the azimuthal field due to the waves. The torsion parameter a is defined by a = (B • V x B)/B 2 , 
and the perturbation of this quantity is 

5a(s,t)=-^-SBt(s,t). (A2) 
RHq 

The vector potential A(r,s,t) is obtained by solving the equations V x A = B and V ■ A = 0. Continuity with the surrounding 
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medium requires that A(R,s,t) is independent of s and t. In the limit An » R, the solution is given approximately by 

A(r,i,0 m y^+-2g- [B*o + ffl *(*.0] ». (A3) 

where r < R and we neglect a small correction A r needed to make V • A exactly zero. The velocity of the waves is assumed to 
be purely azimuthal, \(r,s,t) = (r/R)Sv c j > (s,t)ct), where Sv^(s,t) is the velocity at the cylinder wall. Inserting this expression into 
equation ©, we obtain the electric potential: 

R 2 -r 2 

*(r,j,t)M — B *v^(*,r), (A4) 

where we assumed 77 = and required Q(R,s,t) = to ensure continuity with the surrounding medium. The axial component of 
the helicity flux follows from equation ©: 

H\\ (r, s, t ) = [A x (v x B + V$)] ■s = -^6v 4> (s,t). (A5) 
The average of Tin over the cross-section of the tube is 

H||(*,f) = -^W*>0- (A6) 

Note that the helicity flux is independent of the unperturbed azimuthal field B^q. 

We now compute the helicity flux for damped Alfven waves. In reality the damping would likely involve a turbulent cascade 
of wave energy to smaller spatial scales, but here we use a simple model based on an effective "frictional" damping. The wave 
equations are 

d(S VlP ) = B d(SB^) Syfat) 

dt 4npo Os t 
d(SB^) _ d(Sv^) 

-^r- Bo ^r> (A8) 

where po is a constant density and r is the frictional damping time. The boundary conditions are 6v<j,(0,t) = 5v<p(As,t) = 0, which 
implies d(SB^,)/ds = at the boundaries. We assume the plasma is initially stationary, 6v^,(s,0) = 0. The initial perturbation of 
the torsion parameter, 5a(s,Q), is assumed to be a nearly linear function of position: 



A' 



5a(s,0) = -Act — — —5-^- cos 

k=0 v ' 



As 



2 

where Aa is the perturbation amplitude (we use K = 15). The solution of the above equations is 

AB A _„„ * 4 ./ s \ O 



Aal-^-i), (A9) 



5vj,(s,t)= 7== e 7 ' -j-rsin (nn-£-) — sm(w n t), (A10) 



fc=0 

a: 

(All) 



n 2 7r 2 V As 

£=0 
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cos(cj„f)H sin([x)„f) 



where n = 2k+\; AB^ = BoRAa/2 is the amplitude of the magnetic perturbation; 7 = l/(2r) is the effective damping rate; 
ui n = ^fll~j 2 is the mode frequency; and Q„ = nnvA/As is the frequency for an undamped mode. The above solution is valid 
for 7 < ttva/As. The helicity flux at the mid-point of the tube (s = As/2) follows from equations (IA6t and (lAlOb : 

H\\ (Ai/2,t) = -f(t)R 2 B 2 v A Aa, (A12) 

where 

f(t) = e-v nt^M 2 — sin(wa + i0. (A13) 
^ (2k + Xy-K 1 Lu 2 k+i 

In the absence of damping (7 = 0), the helicity would oscillate back and forth between the two halves of the tube, and the time 
average of Tin would vanish. With damping there is a net transport of helicity from the region with high a to the region with low 

Ti l{ « -f R 2 B 2 v A Aa, (A14) 

where /o is the time average of /(f). We assume that the next reconnection event occurs at a time T which is only a few times the 
Alfven travel time, say, T = 3As/va- Then ft) is given by a time average over the period < t < T. For a damping time t such 
that 0.6 < jAs/va < n, the time average of fit) is given by ft) w 0.01. In section 3 we use the above expressions to estimate the 
hyperdiffusivity in a stochastic field. 
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